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Abstract. The fixed energy scattering matrix is defined on a perturbed stratified medium, and for a class 
of perturbations, its main part is shown to be a Fourier integral operator on the sphere at infinity. This 
is facilitated by developing a refined limiting absorption principle. The symbol of the scattering matrix is 
shown to determine the asymptotics of a large class of perturbations. 

1. Introduction 

In this paper, we study the structure of the scattering matrix on a perturbed stratified medium. In 
particular, we show that its main part is a Fourier integral operator. En route to proving this theorem, 
we develop an improved limiting absorption principle for a large class of perturbations, using techniques 
of Fourier and microlocal analysis. As an application of our results, we prove that the asymptotics of a 
perturbation can be recovered from the scattering matrix at one energy. 

We recall that a stratified medium is a model space in which sounds waves propagate with a variable 
sound speed which depends on only one coordinate. Thus, if wc write the coordinates on K" as z = (x, y) 
with X € R"~^ and y G K, we take the wave speed to be of the form co(j/) and study the wave equation 

where A = - Y.%i ^ We assume that cq is constant for \y\ large and that it is piecewise smooth. We let 

c+ = hm co(y), (2) 

c_ = lim co(y). (3) 

In general, we do not require that c_|_ be equal to c_, but some of our results are stronger when they are 
equal. 

A perturbed stratified medium is a medium on which the variable sound speed, c, has the property that 
c — Co is well-behaved at infinity. Many previous papers have studied the case where the perturbation c — cg is 
rapidly decaying. In particular, precise asymptotics for (c^A— (A— jO)^)~^/, when / G 5(M"), were proved in 

. The inverse scattering problem for exponentially decaying perturbations was studied in |[ |lj, |l^ |9) , 
where it was shown that under certain conditions, the perturbation can be recovered. 

Here we study the case where the perturbation c — Cq has an asymptotic expansion in homogeneous 
terms at infinity. Under certain conditions on c and Cq made more precise in Section |^, we show that the 
scattering matrix for c^A is a Fourier integral operator and describe its singular set. Moreover, we show 
that the asymptotics of the perturbation can be recovered from the scattering matrix at fixed energy. We 
also establish the lead term of the asymptotics for the limiting absorption principle. 

Our results use techniques developed by Joshi and Sa Barreto, |l^, |2l|| to study inverse 

problems in other settings, which build on work by Melrose, [^, and Melrose- Zworski, ]23|, on the structure 
of the scattering matrix on asymptotically Euclidean spaces. As in those inverse results, the fundamental 



idea here is to compute the symbol of the scattering matrix by solving transport equations along geodesies 
on the sphere at infinity. These equations express the propagation of growth at infinity. 

The analysis here is, however, considerably more involved as the unperturbed wave speed cq is not smooth 
on the compactified space achieved by adding the sphere at infinity, even when co (y) is a smooth function of 
y. This is because cq does not have nice asymptotics in |z|. The upshot of this is that cq is well-behaved on the 
compactified space only after the space has been blown-up on the equator at infinity. This manifests itself in 
our analysis by requiring the geodesic flow at infinity to be refracted and reflected by the equator. It was also 
seen in |^ that it makes the asymptotics in the limiting absorption principle much more complicated. There 
is a certain similarity here with many-body scattering, compare, e.g. p5[ |. There the scattering problem is 
complicated by the presence of a potential that does not decay in certain directions and thus appears as a 
spike on the sphere at infinity which causes refractions and reflections of the geodesic flows, Indeed, the 
case where c+ = c_ bears much resemblance to the many-body case. However, when c+ and c_ are different, 
there are effectively different energy levels in the two hemi-spheres, which introduces new complications not 
present in th e ma ny-body setting, and much of this paper is dedicated to coping with those complications. 



In Section O we define the scattering matrix, and its "main part." In case the operator Co(Dy -I- p ) has 
no eigenvalues as an operator on L^(]R, Cq ^dy), then the main part of the scattering matrix is the same as 
the scattering matrix. 

Our first main result is 

Theorem 1.1. Suppose c, co satisfy the general assumptions of Section ^ and either hypothesis (HI) or 
(H2). Then, if 0+ = c_, the main part of the scattering matrix is a zeroth order Fourier integral operator 
associated with broken geodesic flow at time n . If > c+, then the main part of the scattering matrix is a 
sum of Fourier integral operators associated with the mapping 

and the mapping 



1 — c\/c^_ < uo„ 



(a;,w„) 1-^ {-C-Uj/c+,-^l - c^_\uj\'^ / c\) 



if LOn < 0. 

Here, when c+ = c_, the geodesic flow is broken at the equator 0) C §"^^. This can be compared 
to the situation for the Laplacian (||2^), or a perturbation of the Laplacian to an integral power (|^), on a 
manifold with asymptotically Euclidean ends, where the scattering matrix is a zeroth order Fourier integral 
operator associated to geodesic flow at time tt on the the boundary "at infinity." An additional analogy is to 
3-body scattering, where the three-cluster to three-cluster part of the scattering matrix is a sum of Fourier 
integral operators associated to broken geodesic flow at time tt (||2^). Other results on the structure of the 
scattering matrix in n-body scattering may be found in pq| 



Further results on the structure of the scattering matrix are given in Proposition 6.4. 
Our central inverse result is 

Theorem 1.2. Suppose c and cq satisfy the general assumptions of Section ^ as well as either hypothesis 
(HI) or (H2), and n > 3. Then, z/ c+ = c_, the asymptotic expansion at infinity of c — cq is uniquely 
determined by cq and the transmitted singularities of the main part of the scattering matrix at fixed nonzero 
energy. If c+ < c_ , then the asymptotic expansion is uniquely determined by cq and the reflected singularities 
of the main part of the scattering matrix at fixed nonzero energy. 

The reflected singularities are those associated to the mapping (ZU, a;„) s- (— ZU, w„) and, for c+ = c_, the 



transmitted singularities are those associated to the mapping uj i-^ ~uj. Corollary |7T| shows that knowledge 
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of c+, c_, and the singularities of the scattering matrix at fixed nonzero energy determine c, within the class 
we consider, modulo a function vanishing faster than the reciprocal of any polynomial at infinity. 

Following the approach to studying the scattering matrix introduced in psf , in Section |^ we construct 
a parametrix for the Poisson operator. This is a key part of our proofs, as it facilitates an understanding 
of the singularities of the scattering matrix. We work particularly by adapting the techniques of [|l8j which 
are essentially a concretization of the approach introduced in [ p3[ . However, the different behaviour of 
the unperturbed operator CqA in different regions at infinity means that the analysis is considerably more 
involved. 

To pass from a parametrix to the actual Poisson operator, we need a good understanding of the behaviour 
of (A - (A - i0)2c-2)-i/ at infinity, when / e {z)-°°L^{R") and (1 - ^(y))/ G 5(R") for some (j) € Cf'(R). 
In practice, the / for which we apply this will be the error from the parametrix of the Poisson operator. 
When c+ = c_, we can do this by modifying some n-body results of and |Q. However, when c+ < c_ 
these results no longer apply, and we develop new techniques. The essential idea of these techniques is to 
repeatedly develop better approximations with improving smoothness properties. Th us S ection |^ is devoted 



to understanding (A — (A — iO) c ) , allowing us to finish the proof of Theorem LI. In particular, we 
prove the following limiting absorption principle 

Theorem 1.3. Let c and cq satisfy the hypotheses of Section |^ and hypothesis (HI) or (H2). For any 

X € C;?°(§J?^i),, / e {z)-°°l?-i^'^), (1 - 4>{y))f G 5(R") for some G C^{R), we have 

x(z/|z|)(A - c-2(A - tO)T'f = e-^^l^l/l^r("-^)/'ao(z/|z|) + u, 
where oq G C°°(§^-1) and ui G (z)"L2(R") for all e > 0. 

Here C^(S"^^) is the space of smooth functions vanishing in a neighbourhood of the equator and in a 
neighbourhood of {(ZU, G S"^^ : LUn = y^l — c^/ci)}. 

In Section we use a modification of some techniques of ]l9| to prove Theorem |1.2| and apply some 
one-dimensional scattering theory to give further inverse results. 

An announcement of some of these results and an outline of part of the proof can be found in the lecture 
notes [0. 

We are grateful to Fritz Gesztesy for helpful conversations and providing useful references and to Jim 
Ralston for helpful discussions. We thank the London Mathematical Society for supporting this collaborative 
research through its small grants scheme. The first author is grateful to the N.S.F. for partial support. 



2. Assumptions and Notation 
Throughout, z = {x,y) G M"-^ x R. 

Both sound speeds c and cq satisfy < Cm < c, cq < cm < oo. Moreover, co(j/) is piecewise smooth and 
there exists a finite j/m so that co{y) — c± when zLy > yM, with c_ > c+. Moreover, all derivatives of cq are 
bounded except at finitely many values of y. This allows cq to be piecewise constant, for example. 

We require that, away from the hypersurface {y = 0}, c — cq be smooth outside of a compact set K, 
and for simplicity we choose y^ so that K C R"~^ x [—yMjUAi]- Moreover, we make requirements on the 
behaviour of c — cq at infinity. We have, for y ^ 0, 

N 

D^idz) - co(j/)) = d: lA^M'' + o(i^r^-^-i"i) (4) 

for any TV and any multiindex a, where 7j G C^(S"~^ \ {(jjJ, 0)}). Here we use the notation that C^(X) is 
the space of smooth functions on X that have all derivatives bounded. We shall take J at least 2 everywhere, 
although sometimes we shall require it to be larger. Some of our results hold under less restrictive hypotheses. 
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Additionally, we shall often use one of the following hypotheses: 
(HI) J = 2, C-)- = c_, c and cq are smooth. 
(H2) J > 4. 

We warn the reader that the choice of the total space dimension to be n rather than n+1 is in disagreement 
with H] and many other papers on the subject. 

We use the notation (w) = (1 + jwp)"'^/^. Throughout, e shall stand for a small positive quantity and C 
for a positive constant, either of which may change from line to line. 



3. Spectral theory of c§A 

In order to define the (absolute) scattering matrix for c^A, we will need some understanding of the 
generalized eigenfunctions of CqA and c^A, particularly of the space that parameterizes them. Further 
details can be found in, for example, Q ^ 

The operators CqA and c^A are formally self-adjoint on L'^ {W\ Cq'^ dz) and L^(R", c^-^dz), respectively 
and have a unique self-adjoint extension. 

Roughly speaking, the spectral measure of CqA can be given in terms of two kinds of families of functions. 
At fixed energy A, the first is parameterized by almost as for the Laplacian, though the generalized 

eigenfunctions are more complicated. Here §"^^ = {lu — (uj,LUn) G §"~^ : Un ^ 0, aj„ ^ ^1 — c?^/c?_}. 
(Compare |2^, Section 2.1].) 

A second type of generalized eigenfunction comes from eigenvalues of Cq(k^ + D^) on L^(R, Cg^^dj/), if 
there are any. If there are any eigenvalues, let \\{k,) < A^k) < • • • < A^(.^^(k) < c^k^ denote the eigenvalues 
of Cq(k^ + Dy). There may not be any eigenvalues, but if there are, there are only finitely many for fixed 

K and the number grows with k^. Additionally, if k > and Xj > 0, then > 0, as can be seen by an 
integration by parts argument (see, e.g., H, Sect. 2.2]). 

Let K^j be the smallest positive number such that Cq(k^ + Dy) has j eigenvalues for all k > k". Let Kj be 
the inverse of Xj (with the same sign), and let tj ^ lim„|^o A|(k) = c^(k°)^. The {tj} are called thresholds 
of c^A. Let T(A) be the number of thresholds tj less than A^. For < j < T(A), lJ e S"-\ let 

$,(z,A,c.) = e-^(^)--/,(y), 

where fj{y) G L^(R, Cg ^dy) satisfies 

cl{^^^+Dl)f,=Xy, 

and note that (cqA - A^)$j =0. 

At energy level A^, we can parameterize the generalized eigenfunctions of c§A by and r(A) copies 

of The continuous spectrum of c^A is parameterized by the same space as that of c§A. 

Because of the described parametrization of the continuous spectrum at fixed energy, the (absolute) 
scattering matrices of CqA and c^A are operators from L^(S'c) ®i<j<T{\) i^(S"~^) into itself. In js) a 
definition of the scattering matrix is given in terms of the generalized eigenfunctions. Here, however, it will 
be more useful to define the (absolute) scattering matrix using the Poisson operator, which wc shall do in 



Section 4.2 



4. The Poisson operator and the scattering matrix 
The Poisson operator is defined as an operator 

P(A) : C,°°(S^-i) efii^ C°°(§"-2) {z)-y^-'^L\R"). 
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Definition 4.1. //.g = (go, 5i, ■-, 5t(a)) e C^°°(§c) ©Hi^ C°°(S"~2)^ then P{X)g ^ u, (c^A- X^)u = 0, and 
u is determined by its asymptotics at infinity: 



, |_(„_l)/2gjA|z|/c 



30 I o 



T(A) 



Uq 



T(A) 



T/ie functions uq, Uj, u satisfy 



u, 

for any e > 0. Here kj , fj are as in Section I 



u, = |x|-("-2)/2e— ^■(^)l-l6,(x/|a;|)/,(2/), 1 < j < r(A); 



d\z\' 



Proposition 4.2 shows that this expansion uniquely determines the Poisson operator. Definition 4.3, using 



Proposition L3, defines the (absolute) scattering matrix via the Poisson operator. The existence of the 
Poisson operator is proved in Section |6|. 

4.1. The Poisson operator is uniquely determined. In order to show that the Poisson operator above 
is indeed well-defined, we shall need a uniqueness result, for whose proof we shall use Propositio n [l.l| . 

In proving the following proposition, we shall use some results of Weder, H] (See also |l^). We 
recall some of his results below. 

Let A = (— i/4)(z • Vz + • z). We define the commutator [A — A^/c^, A] as a quadratic form (See the 
proof of Theorem 5.4, (28).). By ||2^, Lemma 3.1] for all A > 0, /i > — A^/c'i, there is a compact operator K, 
a compact interval A containing fi, and f3 > such that 

iEA[A- X'^c-'^,A]Ea> f3EA + K (5) 

where Ej\ = E\{A — X^c^^) is the spectral projector for A — X^c~^. 

The following proposition and its proof, included for the convenience of the reader, are essentially adapted 
from |, Lemma 4.17]. 

Proposition 4.1. If u £ (2:)'L^(M") for every e > and (A — A^/c^^)m = 0, then u = 0. 

p8[, there is no nontrivial null space of A — A^/c^, so it suffices to show that 



Proof By the results of ]g7| 

For e,(5 > 0, let us ^ {I + 6{z))-''u e L^{W'). Let L = A - X^/c^, and let $ e C^°°(IR) be 1 in a 



neighbourhood of 0. 
Note that 



where 



us + fc5{z)u 



\fes{z)\<Cj^{l + S{z)y 



(6) 



and the constant C is independent of e and b. 
Then 



[\L,A\^[L)us,^(Ly&) = -2nm(AL$(L)u5,$(i)u5) 



-1 d 



-2jIm(yl<i>(L)(E ^(1 + b{z))-'^us + f,s{z)u), '^{V)us). (7) 



Using this equality, (^), and the fact that 

$(i) : (0)-''i2(K") ^ (z)-''L2(M"), 



we obtain 

\{[L,A]<i>{L)us,<^>{L)us)\<eC^\\u5\\mL)us\\+C2. (8) 

Here and below C$, C2 are constants which may change, independent of e and S, but depending on $, and 
C2 depends on || as weU. Since (1 — ^{L))u = 0, we have 

\\us\\ - mL)us + ML), (1 + S{z))-']u\\ < \mL)us\\ + C2. (9) 

However, using and the fact that is not an eigenvalue of L, we obtain 

{[L, AML)us, nL)us) > MMLhs, HL)us) 

for some /3i > 0, if the support of $ is chosen sufRciently small. 

By choosing e sufficiently small, then, and using (||) and (H), we get that 

PimL)u5f < ylk^f + C2. 

Then 

mL)us\\ < C2, 

and using (0), this shows that for sufficiently small e, \\us\\ is bounded by a constant independent of S, and 
thus u e L^M."), and u ee 0. □ 

We shall use the following notion of an "outgoing" function. 

Definition 4.2. A function u G (z)"'^/^+'i^(R") will be called outgoing if it has a decomposition u = uq + 
Uj + u with the following properties 

u, = |x|-("-2)/2e— 1 < j < r(A); 

for any e > 0. 

Proposition 4.2. Given f e L^(M"), there is at most one outgoing u € (z)-^/^+'^L^(R") with {A— X"^ / Cq)u = 
/. 

Proof. Suppose there are two such u. Then by considering the difference we can reduce this to the case 
/ = 0. Then 




where eo,ej,e, / G (z)'L^(R") for all e > 0. Using the facts that / fj{y)fkiy)'^y = if j 7^ fc and fj is 
exponentially decreasing, this implies that as i? — > cxd 

J\A=R c , J\z\=R , , (11) 

Since Uj £ {y)-°"(xy/'^+'L^{W'') and uq e {yy/^+'(zYL^{«''), we have uqUj G {y)-°°{zy/^+^'L'^{W'). 
Therefore, the right hand side of (pi]), considered as a function of R for large R, is in i?^/^+^'L^(M+), so 
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that uo, Uj e (z)1/^+'L2(K"). This m turn means that 6^ = u-j for 1 < j < T(A). Now we have 
^^it = {-~iX/c)u + u, u £ (z)'^_L^(M"), any e > 0. Repeating the argument almost as above, we obtain 

H'<cf \uu\<c{( HY^'if N^/^ 

l\z\=R J\z\=R J\z\ = R J\z\=R 

and thus 

J\z\=R J\z\=R. 

SO that u G (z)'^L^(R"), any e > 0. By the previous Proposition, u = 0. □ 

4.2. The Absolute Scattering Matrix. In order to use the Poisson operator to define the scattering 
matrix, we shaU need the foUowing proposition, whose proof is a coroUary of our construction of Section ^ 
and Theorem |l.3| (see Section |6^ ). 

Proposition 4.3. Let g = (go, ffi, ffT(A)) G C^(SJ?-i)©^^^^C°°(S"-2). Let9= € K , for K a compact 
set in S"^"'^. Then, for 9 G K, as \z\ oo, 

P{X)g = |zr("^^^/2[e''l^l/^go(0)|K + e-'^l^l/^(ga^))|K] + uk 

where uk e {zYL'^{W^ n {K x [1, oo))) for any e > 0. 

Let y e Ky, Ky C M compact, and let 9 = Then, as \x\ —+ oo, 

T(A) T(A) 

PiX)g = |x|-("-2)/2[^ e-^(^)l-l.g,(0)/,(y)|K^ + ^ e-^-^W\^\g'^{9)f,iy)\j,J+u, 
1 1 
w/iere {t^ € (a;)'L2((M" n ({|x| > 1}) x Ky)) for any e > 0. 

This information about the Poisson operator allows us to define the (absolute) scattering matrix A{X). 

Definition 4.3. The (absolute) scattering matrix A{X) is given, for g e C^(§J?"^) C'^{S"-'^), by 

= ff' G L^i^c^^) ffi ^^(5"""^), where for any compact set K C §c~^ {9o)\k is as in Proposition 
and, for 1 < j < T(X), g'j is as in Proposition [/^. j| . 

We remark that this definition differs slightly from the absolute scattering matrix discussed in 1^. How- 
ever, as the two differ by a straight-forward normalization, we shall use this definition here both to emphasize 
the similarities with the absolute scattering matrix as defined in |22 and because it is most convenient for 
the inverse results. 

For fixed A, A{X) = (yly (A)), < i,j < T{X), with the Aij operators. We cah Aqo{X) the "main part" of 
the scattering matrix. If the operator c^{Dy + p^) has no eigenvalues on L'^ {R, Cq'^ dy) , the "main part" of 
the scattering matrix is just the scattering matrix. 

5. The approximate Poisson operator 

For g — (go, 51, ■••,5t(a)) G L'^(^c~^) ©i"*''^'' i^(§"^^), let Il^g — gj. Since our inverse results involve the 
main part of the scattering matrix, Aoo(A), we are most interested in -Po(A) = P(A)no, where P{X) is the 
Poisson operator. Here we construct an approximation Pq of Pq. 

Let z e M", A e R, A 7^ 0, and let uj e S"^-^, uj = {uJ,ujn), w„ ^ 0. We will construct a (partial) 
approximate Poisson operator Pq{z, X,lu) such that 

(c^A - a2)F e (z)-°°l2(R"), (1 ~ x(2/))(c2A - X'')Pa e 5(R") 

where xiu) G C't?^(II^) is 1 for \y\ < yM + 1, and, so that distributionally as \z\ — > oo. 



Po{z,X,U:) = |zr("-l)/2 |^g.A|.|/c„(,)^^(^)_^g-,A|.|/cofe)^(^^^)^ ^0(1^1 



-(«+1)/2n. 
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We will show how to construct such an approximation when a;„ > 0; the case of a;„ < is quite similar. The 

construction involves solving away errors at infinity. Since the model operator CqA has different behaviour 
depending on the region "at infinity" (y > i/m, \y\ < Vm, or y < — j/m) the techniques involved necessarily 
depend on the region in which z lies. 
We shall need the following 

Lemma 5.1. For f G L'^{[—yM,yM]), the boundary value problem 

cl{X'\w\ycl + Dl)b-XH = f 

with boundary conditions 

-iXuJn/c+b{yM) - b'ijjM) = ai 
iX{l/cl - |a7|2l/4)V26(_y,,) - b'{-yM) = a2 

has a unique solution in L^([0, j/o], Cq ^dy) ifl/c^ > \lJ\'^/c\. 

Proof. This boundary problem can be reduced to the form 

-iXuJn/c+b{yM) - b'{yM) = 

iX{l/cl - p\yciy/H{^yM) ~ b'i-yM) = 

cl{X^\Ljf/cl + D^)b-X^b = g. 

This has a solution if the adjoint operator has no nontrivial null space; and the solution is unique if the 
homogeneous equation has no nontrivial solutions. 
The adjoint operator is the operator 

cl{X'pf/cl + D'^)-X' 

with domain 

{h e i^([0,yo],Co^dy) : -iXuJn/c+h{yM) + h'{yM) = 

and iXil/cl - pf/ciy/^hi-yM) + h'{-yM) = 0}. (12) 
Suppose 5 is a nontrivial element of the null space of the adjoint operator. Then 

/VM 
{cl{X^p\^/cl + Dl) - X^)ggcfdy 
-Dm 

= -g'{yM)g{yM) + g'{-yM)g{-yM) + giymWiyM) - g{-yM)g'{-yM) 

/Vm 
[{X'\lD\ycl-X%')\g\' + \Dyg\']dy (13) 
-VM 

Using the boundary conditions, we find that the second line is equal to 

-2iXcJn/c+\g{yM)\^ - 2iX{l/cl - |a;|V4)'/'|fl(-yM)|'. 

Since the third line of (13) is real, this means that = giyn) = g'iyAi), or g{y) = 0. 

A similar calculation shows that the original operator has no nontrivial null space. □ 

Let u) = {ZUjUn) with w„ > 0. In our construction of the approximation to the Poisson operator Pq, we 
begin with the function ^{z, X, co) which is defined by 

e'^--''/-+<f>+{y) 

where ^+ satisfies 

cliDl + X\l - u>l)cf)^+ = X'^+ (14) 
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and, as y ^ oo, 

<P+{y) ~ e'>^y^-'^+ + i?+(A,c^„)e-*^""^/^+ (15) 

and as y ^ — oo, 



My) ^ r+(A,L.„)e^Wi/-^--i/<+-S/4 (16) 

where when 1/c^ — l/c+ + w^^/c^ < we take the square root so that the right hand side of (^6|) is 
exponentially decreasing. We have (cqA — A^)$ = 0. Note that, up to a constant multiple which depends 
only on n. A, and c±, is the Schwartz kernel of the (partial) Poisson operator Pq.o for CqA when uj is in 
the upper hemisphere of S"~^. We use this as our starting point. 

When y > t/m, we use the techniques of to construct P. Note that when we apply c^A — A^ to $ we 
obtain an error which, for y > j/m, is of the form 

where oi, 02 are classical polyhomogeneous symbols (in \z\) of order —2. We use the techniques of |l8) to find 
a term of the form ae''^^ '^/'^+ with a S 5""^^ which will solve away the first term in the error ( [l7| ) here. Just 

as in [Q, this is done iteratively, solving away an error d G e.^^^^^'^* {S~^gm.odS~j^~^) by using &e*'*'^ "/'^+ 

with b e Spi^^^, and b solving the transport equation 

db 

- 2ic+Xuj ■ Tr = d (18) 
az 

along the geodesies on the unit sphere at infinity. We choose b so that b is smooth at z/|z| = a; in order to 
keep the right coefficient of e*'''l^l/^+ in the distributional asymptotic expansion. 

Let s be the geodesic distance on the sphere from uj and let 9 be angular coordinates about lu. The 
equation ( p^ ) can be solved, modulo an error of lower order, by bi^-j+i\z\~^~^^ just as in ||l^. Section 2], 
giving 

bi^^,+^{s,9,u;)= ! [\s\ns'y-^d_,j{s'Au;)ds' (19) 

zAc+^smsj-' Jq 

where d^jj\z\~^ is the term of homogeneity —j in the error. Note that as long as z/\z\ is in the upper 
hemisphere we are away from — w so the transport equations have a smooth solution. 

We find bi^-j iteratively and then use Borel's lemma to asymptotically sum them, obtaining a bi such 
that 

(c^A - A2)($ + &7e*^^-"/"+) = e^^'^-Wc+g-^Ayc^./c+Q^ ^ 0{\z\-°°) 

when y > yM- Note that the construction of 6/ has not changed 02. 

We will apply almost the same technique to solve away the error e*'^^ '^/'^+e^''^^'^"/^+a2, y > yM, away 
from z/\z\ ~ (— tU, a;„). Here we will use solutions to the transport equation where we choose the initial 



condition at y/\z\ — 0, and the solutions, in analogy to (19), are of the form 



2Ac+(sins)J-i 



(sins')-'' ^d-jM{s' ,0;uj)ds' + Cr^j^i 



(20) 



Here s is the distance on S"^^ from 6 = z/\z\ to the point (ZZJ, — a;„) and 6 is the angular coordinate about 
{u;,—uj„). The value s — sr^ corresponds to On — 0, and Crj depends only on uj and 6. We postpone to 
Section |5.l| discussion of the form of the parametrix near z/\z\ — (— ZU, a;„). 

In the lower hemisphere, we use a similar technique if 1/c^ — |ZUp/c^ > 0. Here the error term is of the 
form 
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where qt € ^phg- Again we have solutions like ( [l9| ) to the transport equation, although this time s measures 



the distance on the sphere from the point (c_a;/c+, ^/l — c'L jcjp/c^). We will have a solution to the transport 

rm 



equation away fr-om {—C-ui/c^, ^ form 



-3 

where 



bT-j+iis,0,uj) 



2Ac_(sins)J-i 



(sins7'~^rf-j,T(s',^;w)ds' + Ctj-1 

To 



The constants (in s) j and Ctj are to be determined. Of course their values affect subsequent errors 
and thus subsequent 6_R,,-j, &T,-j- 

We will use a different technique to construct the solutions when |y| < yi\i. We choose the solutions so 
that Pq is on K". We point out that if cq is not smooth, for example, if it is piecewise constant, we should 
expect it to be impossible to find a smooth Poisson operator on M". 

The values of Cuj , Ctj are determined by solutions to boundary value problems that arise in constructing 
the parametrix when \y\ < hm, as described below. 

When |y| < we look for an approximate solution of the form 



j>0 

The term \z\~^bM,-j is of lower order and is included to improve the regularity at y = iyM- Note that 
(c^A - a2)(|z|-J"6(^, y)e'^"-^/^+) = {cl{Dl + X'\uJ\ycl)b - XH) e^^-^/^+ + 0{\z\-^~^). 



Therefore, for \y\ < yM, to solve away an error of the form \z\ ^ dM.-j{-^, y) we look for bM,-j such that 



X 



4iD'y + X'P\'/c\)bKU-j - X'bM,-j = dMAj^yV)- 

The boundary conditions which bM.j must satisfy come from matching with the solutions in the top and 
bottom hemispheres in order to get a function. They are 

6,^.,(x/|a;|,0)e^^""^-/^+ +e-*^""^-/^+CKj(a;/|a;|)(sins«J-V2iAc+ - 6m,_,(x/|x|, ^m) 

zAc^„/c+[5/,_,(a;/|a:|,0)e^^'^"^-/^+ - e-'^""'^-/=+Cfl,(a;/|a;|)(sinsflJ-V2zAc+] = 6W^_^.(x/|x|, ^m) 

g-.A Vi/ci - PI {x/\x\) (sin sto ) ~' /2* Ac_ = bM.-Ax/\x\, -vm) 

iXa/cl - l/4|cJ|2)i/2e-*^VV^^^^V42'-CT,(x/|x|)(sinsTj-V2*Ac_ = 6^, _^.(x/|x|, -j/m) 

(22) 

where is known (it is determined by integrals over portions of geodesies of {sin sY^^dj) and Crj, 

Cnj are to be determined and are independent of y, and so can be treated as constants in solving the 
boundary value problem. They can be eliminated from this set of equations, resulting in a boundary value 
problem of the type considered in Lemma 5.1, which guarantees us a unique solution to the problem when 
> \uj\/ c\. This then determines bR-j and bT,-j, since Ctj and Crj are determined by bu-j- 
In order to ensure that our function will be at y — yM and aX y = —yu we will add an additional term 
6mj whose total contribution will be of order jzl^-'^^. Let x & C'^(IR), x(t) = 1 for \t\ < 1 and x{t) — for 
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t\ > 2. Let 








= lim I 

yiVM V 






= lim i 

ViVM \ 






= lim 






= lim 


d 
dy 



(e^^^-"/-+(6,,_,(^/|z|) + e-^^^'^"/^+6«,_,(z/|z|))-5M,-,(x/|xUM) 



Note that by om' choice of 6j, Puj, luj, Plj and 7^-,- all have leading order Now, let 

h.-, = (x (^^^^) iPu, + {y- VMhu,] + X (^^^^) IPl, + m,]) H (23) 

For |y| < j/M, this determines the approximate solution of the form (pT|). 

If 1/c?. — |ZZJp/c^ < 0, then we use a slightly different method for finding the approximate solution when 
y l£ Vm- Here, in a manner similar to that used for \y\ < yM above, we solve away the error term by using 
an approximation of the form 

where 

d-j is the coefficient of |2;|~Je*^^ "/'^+ in the error term and is exponentially decreasing in y when j/ < 0, 
and hL,-j is square integrable on (— oo,yAf]. We need in addition a boundary term at ?; = j/m, and this is 
provided by the first two equations of (22). An argument like that of Lemma 5^ shows that there is a unique 
solution to this problem. As in the previous case, &L,-j is chosen to improve the regularity at y = yu- 

Remark . Wc remark that this construction can be carried out, with some minor modifications, for sound 
speeds ci = c + d, where c — Cq has an asymptotic expansion of the type (^, and d is supported in \y\ < yM, 
with d^J2j>a \z\^'dj{x/\x\,y). 

5.1. Approximate Poisson Operator Near its Singularities. For ujn > 0, it remains to describe the 
approximation of the (partial) Poisson operator near z/\z\ = 6 = (— w,a;„) and, if 1 — c^|ZUp/c^ > 0, near 

z/\z\ = 9 = (— c_cJ/c+, ~ c?_ |ti7p/c^). The approximation in these regions contributes to the scattering 
matrix. As these two are quite similar, we will concentrate on the first, using the techniques of jl^. Section 
3]. 

Let w — {wi,...,Wn) — {w',Wn) £ M", and rotate the coordinate system so that uj is the north pole. 
Denote by J^'" the class of operators whose Schwartz kernel can, near the south pole, be written as a 
Schwartz function plus a term of the form 

with a G C^([0, e) x [0, e) x S"~^). From the results of this class is asymptotically complete in 7. 
Moreover, a stationary phase computation which can be found in jist shows that away from the south pole 
this is equivalent to the class of operators whose kernel is of the form e''^'""6, with b a polyhomogeneous 
symbol in |w| (and smooth in w/jwl) of order —7 + (n — l)/2. We recall below some additional facts about 
the operators I^'" from [ p8[ . 
We recall from 

11 



Proposition 5.1. Ifu{w,uj) G I^^" and f € C°°(§"-i x §"-i), then 

e*^l"'l y u(|u>|6i,w)/(6i,w)d6idcj 

is a smooth symbolic function in \w\ of order —1 — a and its lead coefficient is \w\^^^"{K, f) , where K is 
the pull-back of the Schwartz kernel of a pseudo-differential operator of order a — 7 — (n — 2) by the map 
9 i-^ —6. The principal symbol of K determines and is determined by the lead term of the symbol, a{t,S,fi) 
of u as S ^ + . 

/^From [|8[ Propositions 3.1 and 3.2], we liave 

Proposition 5.2. If u E I^'" and V ^ J2j>2 then 

(A - T^)u e 



^From Lemma 3.2 of [ |18[ 

;-7,("-3)/2 



Lemma 5.2. If u e /^'^"-^^/^^ y ~ J2j>2 H^^Vj{w/\w\), and (A + V - t^)u e iV ' , then (A + F - 
Again from 

Proposition 5.3. If u £ /^'" = n^/^'", then u — e'"^^^ f{w), with f a classical symbol of order —a — 1. 



We proceed as described in 1 18, Section 3]. Using the first part of the construction, we have an approxima- 
tion of the Poisson operator that blows up as — z/\z\ approaches (— ZU, ujn) or (— c_cJ/c+, — y^l — c^|iUp/c^). 
Recalling that lo^ > 0, near 9 — (— ZU, a;„), we find an approximation of the Poisson operator of the form 
u7?.+ , where u e -^a/c+ ^"^^^'^^ ^y^^ g^j^j . (ij^ ~^\)- For 9 near (— c_ZU/c+, — y'l — c^|ZUp/c^), 

the approximation of the Poisson operator is of the form uTZ^, where u G -^a/" '^''^^ and 7?._ : 



(a;,a;\) 1-^ (J_cj/J+, y 00 — Ji|w|^/J^). Putting aU of this together, we get an approximation Pg to the 

(partial) Poisson operator with a remainder term (c^A — X^)Po that is in {z)~°°L'^, and is Schwartz after 
multiplication by a function </> S C^(M) which vanishes for \y\ < i/m + 1- This error is then solved away by 
applying (A — (A — iO)^c~^)~^c~^ , and it is for this that we need the results of Section |[ Taken together, 
the approximation of the (partial) Poisson operator described in this section, together with the results of 
Proposition and Theorem 1.3, show Theorem |l.l| . 

5.2. Approximation of P{X)Ilj, I < j < T{X). For completeness, we briefly outline how to construct an 
approximation Pj{X) to P{X)Ilj, I < j < T{X). The approximation will have the properties 

(c^A - X^)P^ e (z)-°°L\K"), (1 - (t>iy))ic^A - A^)^,- e 5(R") 

if G C^(R), (f>{y) = 1 when \y\ < j/m + 1, and, distributionally as \z\ —f 00, 



P,{z,X,ZJ) = |xr("-2)/2 + 



-0(|x|-"/2(y) — ). 



Here uJ e S"-^ 

For the construction of Pj, we begin with $j = e^'^'^^'^^''^ fj{y), which is, up to a constant multiple 
depending on n and Kj, the (jth partial) Poisson operator for CqA. We have 

(c^A - A^)*, ~ kr'e-^(^)--d^fc(:E/|x|,y), 

fc>2 
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where d^k is smooth in exponentially decreasing in y, and smooth in y when \y\ > yM- 

To solve away the error with k — 2, we write 

d-2{x/\x\,y) = d_2,i{x/\x\,y) + d^2,2{x/\x\)fj{y) 

where / d^2,i{x/\x\,y)fj{y)cQ'^{y)dy = 0. Since (i-2,1 is orthogonal to fj, we can find g-2 such that 

(cliDl + nj) - \')g_2ix/\x\,y) = d^2Ax/\xly) 

and use e'^'^^^^'>'^''^\x\~'^ g-^2{x /\x\, y) to solve away the e^'^^'''^^^'^\x\^^d-2,iix/\x\,y) error, up to a term van- 
ishing one order faster at infinity. We note that g^2 is exponentially decreasing in y, since c?_2,i is, and that 
this term does not contribute anything to the scattering matrix, since distributionally it is 0(|a;|~(""'"^^/^). 

To solve away the error \x\~'^e^'^^^'^'>^'^d^2.2{x/\x\)fj{y), we use the techniques of in the x variables 
only. That is, essentially as in our construction of Pq, we solve transport equations alon g ge odesies on §"^^ 
beginning at x/\x\ = to. Near x/\x\ = —uj, we must use the second ansatz, as in Section 



5.1 



or 1^ 



The subsequent errors are solved away in exactly the same manner, resulting in an approximation as 
claimed. 

6. End of proof of structure results 
In order to finish the proof of Theorem |1.1| , we must tie up a number of loose ends. These include showing 



the existence of a function u as in the definition of the Poisson operator (Definition 4.1), proving Proposition 



4.3, and showing that the construction of the (partial) Poisson operators in the previous sections captures all 
the singularities of the scattering matrix. As these use similar techniques, we give the proofs in this section. 
For A € M \ {0}, e > 0, the limit 

iim(z)-v2-.(,2^ _ _ z5)2)-ic2(z)-v2-. = iini(,)-i/2-(A _ - ^5fr^zry^-^ 

oJ.U oil) 

as an operator on L^(R") exists in the norm topology (|^, |^). In this section, we study further properties 
of (A - (A - iO)2c-2)-i when it is applied to a function / G {z)-°°L'^{W), (1 - 4){y))f G 5(M") for some 
4> G C^(R). In particular, we are interested in the asymptotics at infinity of the resulting function. 

For simplicity of exposition, we shall assume A > throughout this section. The results for A < can be 
proved in a similar way. 



6.1. "Outgoing" solutions. The main result of this section is Proposition 6.1. This proposition, combined 



with the approximation of the Poisson operator in the previous section, shows the existence of a function u 



with the properties given in Definition 4.1. Taken together with Proposition 4.2, this shows the existence 
and uniqueness of the Poisson operator. 

In proving this proposition, as well as in many others, we shall use the fact that if 

u = lim(A - c-2(A - ief y^f - (A - c-'^{\ - iQf T^f, 

then 

u = {^-{\-iQfc^^)'\Vu + f) (25) 

where 

V^\^{c~^~c^\ (26) 

Additionally, 

(A - (A - ^Qfcf)-^g{z) = (2^)1- j e-«((i?^ + - (A - iGfc^^y'^. ■)){y)di (27) 
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where g{£,,y) = J e~^^'^g{x,y)dx is the Fourier transform in the x variables only. We shall repeatedly use 
this notation for the Fourier transform in the x variables only, and V is given by (|2^). 

We make several remarks about the operator {Dy + — {X — iO)^CQ^)~^ . As an operator from {Dy — 
(A - iO)^c^^)-^ : (y)-i/2-<^L2(M) ^ {yy/'^+'L'^{'R) for any e > it is smooth for \t\ < A/c+, t away from 
A/c±. Near t = A/c+, (AV4 - t2)i/2(£)2 _,_ ^2 _ (\ _ iQ)2^~2yi ^ smooth function of (A^/^ - i2)i/2^ 

and, if c_ > c+, near t = A/c_ it is a smooth function of (A^/c^ — t^)^/^. 



We shall use the following lemma in the proof of Proposition S.l 



Lemma 6.1. If u — limej^ol^ — (A — ie)'^c and f E (z) °"L^, then Vu{^,y) and all its derivatives 

with respect to ^/|C| are, for e > 0, m "'-1/2-^-/3 (R^-i. {y)~l3 L'^{^y)) for < f3 < J - 1/2 - e. 

Proof. Throughout the proof e > is small, and may change from line to line. 
We have (A - X^/cl)u = + /. Then 



Vu{i,y)^C e~^'^<V{x,y) / e^^-^{{Dl + \^\^ - cl{X~ zQf)-'[f{i^,-) +Vu{r^,-)]){y)dvdx. 



(28) 



Let ZJ^j/i^Dj stand for a derivative tangent to |^| =constant. Then 

{D^aii^^y J e— «y(x,y) I e^-^{{Dl + M'-c^\X-iOf)-'Vu{v,-)){y)dr^dx 



+ E 7 y, j e-^Cp^ + |^|2 „ ^(A - zQfy^ V^^h ■)){y)dvdx. (29) 

\a\<j 

Here D'^ri/\ri\) ^ derivative of order \a\ in the variables and 

\D^,D>;DjV^{x,y,0\<C{l + \z\)-''\^\~'^ (30) 

for large \z\, where if fc > 0, we need also require \y\ > yM- That is, Va is a symbol in x, a property we shall 
use. Let A^/|^| be the Laplacian in the ^/|^| variables. If < t < 1/2, then 

(A^/I^l)* / e-^-<V{x,y) j e^--\{Dl + - c^^X^)-^Vu{r^,-)){y)di^dx 

e'^'^<V{x,y) j e^^-^iDl + H-" - c^^X^)-^ Vu{r^,-)){y)di^dx 

+ j e— «T4(x,y,e) j e^'^-\{Dl + \r^\'~c^^\'r^Vu(T^,-)){y)di^dx (31) 

where Vt{x,y,£) satisfies, for large \z\, 

\D^D^^DlVt{x,yM < C„,t,^,„(l + \,\)-J-i+2t-\a\-m 

where we require m = if |?/| < t/Af . 

Because of the decay properties oiV, Vu{v,y) e H-^-^/'^-''-^{Wy^; {y)~^ L'^{Ry)) for J-e-1/2 > /3 > 0. 
Of course, /(ry,t/) e i/-(R;;-i; (y)-°°4). 

Suppose for some t > 0, {A^/i^iY/^Vu G i7^-i/2-£-/3(j^n-i. (^y)-P L^(^Ry)) for J - e - 1/2 > /3 > 0. Then 
using equations (H), (29) and (31) we obtain (Aj/|^| )(*+2^-2-2.)/2^ ^ iZ -'-1/2-^-/3 (r"-!. (y)-/3L2(Mj^)). 
Since J > 1, the proof follows by induction. □ 
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Lemma 6.2. Ifw,!]^ M."\ and 

\D:D^^aiw,rj)\<C^f,{l + \w\r\^\ 
for all multiindices a, (3, then the operator A defined by 

(A<?)(7?)= J ^^■^aiw,ri)g{w)dw 

is a continuous map 
Proof. We can write 

{Ag){r])^ I e''^-^'^-'^'^a{w,ri)g{ri')dri'dw 



with g G H '(R™) the inverse Fourier transform on R™. Since a is a symbol of order in w, this is 
a pseudodifferential operator acting on g and the result follows from standard pseudodifferential operator 
theory. □ 

Lemma 6.3. Let a(w,ri) be a symbol of order in w = {w,Wm+i) G R™"*"^, and let h{g) E _ff~'''(R'") be 
supported away from t^ = \g\'^ , where t is a nonzero constant. Then 



provided 7 > 0. 

Proof. Let g e {wy^^-^/^ L^{R"^+^). Then 



(g, Ah)^iJ aiw, ry)e-™"'-™'"+^-g(z«)du;|^^^p3^, /i(ry)). (32) 
By the previous lemma and the restriction properties of elements of Sobolev spaces, 

J a(«;,,y)e-™"'-™-+^-g(z«)du;|^^^^5-^ € H'< 
if 7 > 0. The pairing (|3|) is then well-defined for aU such g, and Ah e {wy+^/'^L'^{W^+^). □ 



Proposition 6.1. If f £ (z) °^L^{W" ) and J >2, then u = (A — (A — iO)^c ^) ^/ is outgoing in the sense 
of Definition {.i. 



Proof. We use (|2^) and \f^. Choose * e C;?°(R) to be 1 for |^| < A-^/c^ and supported in a slightly larger 
neighborhood. We use the fact that we can write, for 5 > 0, (1 — ^{Ol){{D\ + — (A — J(5)^Cq ^)~^ as a 
sum of an operator bounded on L^(R) and an operator involving projection onto the discrete spectrum. The 
part corresponding to an operator bounded on I? gives an element of I? . For the part corresponding to the 
discrete spectrum, near the poles we move the contour of integration in |^| slightly into the upper half-plane 
when a; • ^ > and into the lower half-plane for x • f < 0, obtaining 

T(A) 

e-«(l - ^{m{{Dl + - (A - t^fcr^^)-\Vu + /)(e, ■)){y)di ^ ^ ^■(^)l-l6,(^)/,(y) + 

where ui, g^ui € £^(R"). We remark here that bj e C°°(S"^^), thanks to the fact that it originates from 
a stationary phase applied to functions smooth in (Lemma |6.l[ ) at |^| — Kj(A). We refer the reader to 
the proof of Theorem 4.1 of Q for greater detail in a very similar computation. 
When \y\ < C for some constant C, we have that away from |^| — A/c+, 

y) = "^{miDl + leP - (A - ^Qf)-\Vu + m, .))(y) e L^i^r') 
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with a norm independent of y, \y\ < C . Near |^| ~ A/c+, y) G LP{W^ ^) for any p < 2, again with norm 
independent of y, \y\ < C, so using the mapping properties of the Fourier transform we have 



When y > yM, for |,^| < A/c+ we have 



{{Dl + \e - (A - tQ)%')-'Vu + c^fi^, ■)){y) = h,{y)e-^V^^7^^ + h^iy, 0- 

Since h2{y,S,) £ L^(R'^^^; (y)^ '^^'^/^^'^ (Ry)) , taking its inverse Fourier transform in a; gives an element of 
(?/)^-'+^/^~'^L^(R" n {y > yAi}), and similarly if we take a radial derivative. 

For the term /ii(y)e"*V^^/^+^^, note that away from |^| = A/c±, hi G i/"^"i"'=(M""i). Notice that 



+ zA/c+)e"-«-^^V^Vc^+-l«P = • e - i^y^AVc^ - |^|2 + zA/c+)e"-«-*^V^Vc^+-l?P^ 



(33) 



that is, it vanishes to first order on the critical points of the phase. Therefore, smoothly restricting our 
region of integration to |^| < X/c^ — 5 and away from |^| = A/c_, an integration by parts argument shows 
us that we get an element of (z)^L^(IR" fl {y > ?/m}) for any e > 0. 

Near |C| = A/c+, \^\ < A/c+, = (A^/c^ - m-'/'b{^ X^cl - mh+iO, where b e and h+{0 = 



Y A^/c^ — l^p), G _ff '■'^-'^/^"'^(M"). Using a partition of unity, we can work on coordinate patches 
on which ^ for some i. For example, let ^ = (^2, ^2, ••■Cn-i) and for > 5 > 0, we use (^,< = 
A^/c"^ — l^p) as our variables of integration. Then we need to compute 

{^+^X/c+) I e^(-V^V4-t^H«l=+^-«-yt)^^(^^i)5(^)^^^( A2/,^_i2_|J|2^J^i)rfi. 

Here xi is a smooth cut-off function and x = (2:2, 2:3, a;„_i). The resulting integrand vanishes at the 
stationary points of the phase, so that we may integrate by parts. Applying Lemma 3.3 and a related 
result for the boundary term, we obtain an element of (z)'^L^(R"). We can repeat this procedure, using the 
partition of unity, to cover the integration over the region with |^| < A/c+, |^| near A/c+. 
For 1^1 > A/c+, we have 



l«l>Vc+ 



a{^W-X^/c\){Vu + f){U)<t>{v'. - \Vc\)dy' + h^^it y)) (34) 



with a{t), (j){y' ,t) smooth in t and \(j){y' ,t)\ < Ce^ *. The function (j){y, y |CP ~ A^/c^) is in the null space of 
Dl + I^P - X^/cl Since /i26, Dyh2b G H-'-^-"-^ {W'-^; (yy^L^iR)) for J - 1 - e - 7 > 0, its contribution 
is in {zyL'^{W^) and can be neglected. As in the previous integral, for the rest we use a partition of unity 



to work in sets with 7^ 0, some i. Near ^1 7^ 0, we use coordinates £,,t = y — A^/c^ almost as before. 

Using the technique of the proof of Lemma 3.3, applying (^^^ + iA/c+) when y > yj^i, integrating by parts, 
and using Lemma |6.2|, we need to check only that 

Wi{Lt,y) r e-y'(f>iy\t)w2{lt,y')dy'dtdl (35) 

t>OJyM JyM 
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is finite when wi G i/'=/2(M"-i; {y)-'^-'^/^L^(R)) and W2 £ H-'-^-^^ {W'-'^; {y)-^L^{R)) for any S > 0, where 
W2 is supported near t = and away from |^| = A/c+. Since by Holder's inequahty 
ry 



e 



' Vm 

for any 6 > 0, using Holder's inequality again we see that (|3^) is finite. 

Near |^| = A/c_, hi = h^i + h_2, where /i_2 G i?'''^^"'^ and can be treated by an integration by parts as 



before. However, h_i is the restriction of an element of H'^ ^ to the hypersurface y A^/c'i — 

Here we introduce the variable of integration t = A^/c^ — this results in an extra factor of t in the 
integrand and makes t = a stationary point in of the phase. However, the extra factor of t in the integrand 
is enough to allow us to integrate by parts after applying + iA/c+, and arguing as in Lemma 6.3, we 

obtain an element of (z)'^L^(R" n {y > um})- 

A similar technique shows that (g^ + iA/c_)(A - (A - iO)^c-^y^ f^y^^y^^ G {zyL'^{W' n {y < -Vm}) 
for any e > 0. □ 

For future reference, we remark that the proof above has shown 

Corollary 6.1. If f e (z)-°°(R"), then for yeK,KcR compact, \x\ > 1, 

T(A) 

(A - c"2(A - *0)2)-VW||.|>i.,eK = |xr("-2)/2 J2 e-'''''^'^^%ix/\x\)f,{y) + m 



1 



where bj £ C°°{§"~^) and ui e {xyL^{W'-^ x K). 



6.2. Smoothness of leading order coefficient. To finish the proof of Theorem 1.1, we need to prove 
Theorem 1.3, which we recall: 

Theorem . Let c and cq satisfy the hypotheses of Section |^ and hypothesis (HI) or (H2). For any x G 

C^iS'^-^), f e {zy°^L^{W"), (1 - </>(y))/ e 5(M") for some (j) G C^{R), we have 

X{zl\z\){^ - c-2(A - *0)2)-i/ = e-'^l-l/-|zr("-i)/2ao(z/|z|) + wi 
where ao € C°°(§J?-i) a«d ui e (z)^^^^]^™-) ^ > g. 

The proof of this theorem will occupy the next two subsections. Before we prove it, however, we give 
some applications. 



Proof of Proposition ^ We have P(A)Hj = Pj{\) - (A - c-^{X - iO)^)-'^{A ~ X'^c-^)Pj{\), where Pj are 
the approximate (partial) Poisson operators constructed in Section ||. Proposition i.2 now follows from our 
construction of Pj , Proposition |6.l| , Corollary and Theorem |l.3| . □ 

Moreover, if we take / = (A — A^c~^)Po(A, z,uj), where Pq is the approximation of the (partial) Poisson 
operator we have constructed, then this theorem shows that the the only singularities of the main part of the 
scattering matrix come from the structure of the approximation of the Poisson operator we have constructed. 
This then proves Theorem 1.1. 

in one case. 



The following proposition implies Theorem 1.3 



Proposition 6.2. Suppose c and cq satisfy the hypothesis (HI). Let x ^ C'^{^'): ^) and let f E (z) 
(1 - 0(2/))/ e 5(R") for some cj) € C^{M.). Then 

N 

x(z/|z|)(A - c-2(A - iOf)-^f = |^|-("-i)/2e-»^l-IAo(y) ^ \z\~^ a,{z /\z\) + 0{\z\-^'''-^'^'^-^-^) 





where a, e C^^^W^^). 
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Proof. Because c+ = c_ , we can recast this in the form of a perturbation of a simple n-body problem and 
use the fact that much is known about asymptotic expansions of the resolvent applied to Schwartz function. 
First, we note that 

(A - c-2(A - tOf)-^ = (A - A2(c-2 - c;2) _ c;2(A - iO)^'. 

The operator A — X^(cq^ — c^^) is a particularly simple example of a class of n-body operators widely 
studied. The operator A — A^(c^^ — c^^), while not quite an n-body operator since the potential depends 
on all variables, is a perturbation that has many of the same properties we desire. 
The paper |^^, which builds on results of shows that 

x(./|z|)(A - A2(co 2 - cf) cfiX - ^Or)-'f 

N 

= |;^|-("-i)/2e-A|.|/co(,) J2 \z\'^b,{z/\z\) + 



with bj smooth. The proof is such that the results of ||24| ho ld with cq replaced by a sound speed c of the 
type considered here. Roughly speaking, this is because [13| requires that the operator A + Vi (where for 
us Vi = — A^(cq ^ — c^^) or Vi = — A^(c^^ — c^^) ) satisfy a Mourre estimate and some regularity and decay 
properties, both of which are satisfied for either Vi. Vasy remarks already that the results of Section 2], 
which are local versions of results of |Q , will hold in our case. Then the results of J^, Section 3] hold for 
our case, since just as in that paper we can argue that from the results of and [[22| that 



and the asymptotic expansion follows from Proposition 2.8] and the remarks made there. (Roughly 
speaking, WFsc(u) provides a microfocal description of the lack of decay of u- see |Q for a definition of the 
scattering wave front set, WFsc, and ■) □ 

Before giving the proof in the case c_, we give some explanation as to why the proof of [2^ does 

not apply in this case. The argument of |24| uses very strongly the fact that, for f £ S, 

where x S Cf>(§J?-i). 

This is not, however, true in general when c+ 7^ c_. The results of Theorem 4.1] show that in general 
for / e 5(R"), 

((A - (A - ^OfcoY'f) (^) - e-*^l^l/^+|zr("-i)/2ao(z/|z|) 

+ |zr("+i)/' (e-'^l"l/"+ai(z/|z|) + e-*^l^l/=— + 0(|z|-("+3)/2) ^gg^ 

when < e < y/\z\ < (1 — c^/c^)^/^ ~ e and \z\ 00. If hi 7^ 0, then the scattering wave front set is not 
contained in R1^/^_^- This means that the scattering wave front set of x(A — (A — «0)^/c^)~^/ is in general 
more complicated than the c+ = c_ case (and unknown, to the best of our knowledge) , and the techniques 



of |2^, 24 cannot be immediately applied. Similar differences can be seen in the resolvent estimates of |1J 
Theorem 1.1] for the n-body problem, and [ p^ Theorem 3.1], for a particular stratified medium. 
Instead, we take a different approach. 



6.3. Proof of Theorem 1.3 in case hypothesis (H2) holds. In order to prove Theorem 1.3 when 
hypothesis (H2) holds, we will make heavy use of equations (^5|) and ( p7[ ) . We use the fact that the more 
rapidly g decays at infinity, the more we can say about (A — (A — iOpCQ^)^^,g, using equation (p7|). To 
take advantage of this, roughly speaking, we find approximations w of u — (A — (A — iO)^c~^)~^ f so that 
(A — X'^/cq)w — Vu + e, with the error e decaying faster than Vu does, and w decaying faster than u. Then 
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u = w — (A — (A — iO)^CQ ^)^^(/ + e) (compare (p5|)). The better rate of decay of e improves our knowledge 
of u. 

In practice, the proof is somewhat more complicated. We study the behaviour at infinity of x(^/k|)u, 
X G C^(§"~^), and we introduce a "microlocal" cut-off '^{Dx) so that X^{I^x)u has the same leading 



behaviour as xu at infinity (Lemma 3.4), but is easier to understand. 

Lemma ||4| shows that, for suitable -^{l - ^[Dx))u G (z)^~^/^L^(R"), so that x ^{Dx)u captures 
the leading behaviour of xu. Lemmas 6.5- 3. £ are preliminaries, and Lemmas 6.10| - ^.li] are used to make 



successive approximations of xw in the proof of Proposition p . 3| , which shows that x^(^x)'^ has an asymptotic 
expansion with smooth coefficient in the leading order term. 

The next lemma shows us that, for suitable 4' and x? X^{Dx)u is the leading order term of x."^- 

Lemma 6.4. Let x e C;?°(§^"i), e C^(M""i), and suppose that if {x/\z\,y/\z\) e suppx and y > 0, 
then ±Aa;/(c+|z|) ^ supp(l — ^), and if {x/\z\,y/\z\) G suppx and y < 0, then ±Ax/(c_ |z|) ^ supp(l — ^'). 
IfJ>A,u= {A^ {X-iOfc-^)-\f, and f e {z)-°° L^iW), then 

xiz/\z\)il-^{Dx))u, -^^xiz/mi-^{Dx))ue{z)-'/'+^L^R-) 

for any e > 0. 

Proof. We give the proof for x supported in y > 0, as the proof for x with support in ?/ < is quite similar. 



The proof closely resembles that of Proposition 6.1. Let 'J'l e C^(]R) be such that supp5'i(t) C {\t\ > 



A/c+}, supp(l — ^fi) C {\t\ < A/c+ + (5}, some small S > 0. Then, by the same type of arguments as in the 
proof of Proposition |6.l| , since the eigcnfunctions of Dy + Cq^A^ arc exponentially decreasing in y, we have 

Xiz/\z\) I e-«(l - ^{0)^,imDl + Id - Co 2(A - zO f)-\Vu{t •) + f{t ■)){y)di e {z)-'/'+^L\R-). 
Moreover, for y > yM, using (c3), we see that when z is restricted to the support of x(-z/|-z|) there are no 



stationary points of the associated phase x- — y^ A^/c^ — on the support of 1 — VE'(^). Therefore, when 
the integrand is supported in \^\ < \/c^ and away from \S\ = A/c±, we can integrate by parts as in the 



proof of Proposition 6.1 and get an element of (z) ^"+^L--(M"). We need only examine the integration near 



Id = A/c± more closely. 

If c-|_ < c_, near ji^j — A/c_, y > yM, we have 



{Dy + \e - (A - lOfr'iVui^, •) + /(d •))(y) - e-'yV^^7^^{a,iO + {\'/cl - l^Y^'a^iO) + b{y, 

(37) 

where ai,a2 E i7'(M""i), / an integer with ? < J - 1, and 6 e L'^{R"-^; {y)-''+^+'L'^{M.)). Since again 
there are no stationary points of the phase x ■ ^ — y^^J A^/c^ — with z/\z\ in the support of x a-nd ^ in 
the support of 1 — 5*, we integrate by parts. Since, for any G C^(R), the Fourier transform in x of 
^ad^l - A/c+)(AVc?. - \£.\^)-^/^ is in LP(R"-1) n C°^{W'~^) for any p > 2, we have the contribution of the 
integration of Ui terms over this region is in (z)~^/'^+^_L-^(R") for any e > 0. Clearly, the contribution of the 
integration of the b from (||) is in (z)--^+^+''L'^{R"). 

Near |d ~ A/c+, as in the proof of Proposition |6.l| , we divide the integration into two pieces. When 

Id < A/c+, we introduce the coordinate t — ^ A^/c^ — as in the proof of Proposition 6T and can 

integrate by parts to get an element of (z)~-'^/^+'^L^(R"). 

For the part with |d > A/c+, we use (34). Changing coordinates as in the proof of Proposition ^j] , 

we can, in the region with > (5 > 0, consider the phase to be xi \Jt'^ + ^ A^/c^ +x-S^+ iyt, using 
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t = Y l^p — A^/c^ and the notation of the proof of Proposition 6J. Again we have no stationary points and 

integrating by parts gives an element of (z)~-'^/^+'^i^(R"). 

d 

d\z\ 



A similar argument gives the same result for gf-Tx(l — '^{Dx))u. □ 



We use the notation S![. ^ = {a; = {uj,ujn) G §" ^ : ±w„ > 0} and will use differential operators of the 
following type. 

Definition 6.1. We say that a differential operator P G DiffJ,(]R") if P is a differential operator of the form 

± \a\+j<l II II II 

where bj £ C°°([0,oo)), aj,a.,± € Cr(Sr^) D'^/\z\ 

is a differential operator of order \a\ in the z/\z\ 

variables. 

We shall make use of the following lemma, whose proof follows by a straightforward computation. 

Lemma 6.5. If P e DiffJ,(M") and b e C°°(§"-i), then [g|y,P] G Diff[,(R"), [g^^,P] G Difft(M"), 
[|z|-"6(z/|0|),P] G |z|-'"Diff^-i(R"), and [A, P] G Diff'+i(R"). 

In order to construct the desired approximation, we shall use the following lemma. In practice, when we 
apply this lemma, the first term will be used in solving away the error, and the subsequent terms will be of 
lower order. In particular, in applications h will vanish faster than g at infinity and so will x±(apj + 7^)9- 

Lemma 6.6. Suppose (A - X^/cDg = h, P e Diff^(R"), and x±,v G C°°(§"-i) with the support of x± 
contained in S^^^. Then 

(A - X^cl) {\zr+'vx±{z/\z\)Pg) = {2j - 2)\z\-^vx±—Pg + 

c± 

(2j - 2)\z\-^vx±P{j^^ + —)g + \z\-^-\x±{Pi+i9) + \zr^'yo{vx±) ■ VoPg 
o\z\ c± 

+ \z\-'+hx±Ph + A{\z\-^+\x±)P9 

= i2j -2)\z\-^vx±—Pg +\z\-^-'Pl+^g + \z\-^+'vX±Ph (38) 
c± 

with P;+i,P/_^^ G DiffJ,+\M"). Here Vq stands for the gradient on S"^^ 

Proof. The proof follows from a straight-forward computation, using Lemma |6.5| . □ 
If G C°°(M"-i), we use the notation 

{^{D,)f){z) = (2^)1-" y e'^'-'--'yi^{Of{x',y)dx'd^. 
Note that [^{0^:),^] = and co(y)] = 0. 

Lemma 6.7. If xi, X2 e C^{S"~'^) with supp xi n supp X2 = 0, * G 5(K""i), and 

Au ^ Xiiz/\z\)'i'iD,)x2iz/\z\)u, 
then A : (z)°'L^{W) {z)-°"L^{W') for any a G R. 

Proof. The crucial observation is that if (x, (a;, y) I G suppxi sjnd {x' ,y)/\{x' ,y)\ G supp X2, then > 
(3\{x,y)\, \x - x'\ > /3\{x',y)\ for some (3 > 0. 
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The Schwartz kernel of -^{D^) is given by (27r)i-"*(a;' - x), where # G 5(M"~i). Thus, for / e L2(M"), 
\Af\{z)^{2ny~^\Xi{z/\z\) J ^{x'-x)x2{{x',y)/\{x',y)\)f{x',y)dx'\ 

< C\xi{z/\z\) Jil + \x- a;'r™-"(l + \x - a;'|r+"#(a;' - x)x2{{x' , y)/\{x', y)\)f{x', y)dx'\ 

< C{zr"\ J{l + \x- - x)\'dx'y^\ J {{x' ,y)r'^\fix' ,y)\^dx')'^^ 
<C{zr"\j {{x' ,y)r^"\f{x' ,y)\^dx'f'^ 

for any m (where the constant depends on m and a) and thus it follows that Af E (z)^'^L^(R"). □ 



Lemma 6.8. J/^- e C^{R"-^) and g E (z)-^L2(M"), then'^{D^)g g {z)-''^L^(W'-). Suppose D^g, D'^Pg g 
(z>-'3i2(K«) p g Difr^(R") and for all multiindices a. Then D<^P-^{D^)g g (z)-f^L^{R'') for all 

P g DiffJ:(M") and all multiindices a. 

Proof. To show that if 5 g {z)- 1^ L'^ [W^) , so is '^{Dx)g, we take the Fourier transform of '^{Dx)g: 

J^{^{D,)g){r^,T)^^{f^)T{g){v,T) 

where here we are using !F{h){ri,T) to denote the Fourier transform of h in all variables. Then, since 
T{g) g i7^(M"), T{^{D^)g) g and 'l![D^)g g {z)-^L'^{W'). 

We give an indication of the proof of the remainder of the lemma. Suppose -D"^, D^Pg g (z)^'^L^(R") 
for aU P g Diff.^(R") and for aU muhiindices a, and fc > 1. For x+ e consider 

(27r)"-ix+(^/|^|) (^|z|A + !M^$(z?,)5 = x+(z/|z|) j i>{x' -x){x' ^,+y^^+^\\{x' ,y)\/c+)f{x' ,y)dx' 
+ i\/c+x+{z/\z\) I ^(x' -xmx,y)\-\{x',y)\)f{x'.y)dx' 



X+{zl\z\) j Y^D^ix' -x)^fix',y)dx' 



Let X e C^{^+^^) be 1 on the support of x- Then, using the first part of the lemma, 

X+{zl\z\) j *(x'-x)x((x^y)/|(x^y)|)(x^V.-+^/|- + ^A|(x^y)|/c+)/(a:^y)d:rg (z)-^L2(R") 
since {x' ■ Vx' +2/^ + y)|/c+)/(a;', y) g (z)^'^L^(M"). By the same reasoning, 

X+{z/\z\) j Y^D^{x' ~x)-^f{x',y)dx' g {zyf^L^W^). 



Moreover, by Lemma 6.7 



X+{zl\z\) / ^(x'-a;)(l-x((x',y)/|(a;',2/)|))(a;'-V,, +y^ + zA|(x',y)|/c+)/(a;',y)dcc 

g (z)-°=x2(M"). 
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To finish, note that {x — x')'^^{x' — x){\{x,y)\ — \{x' ,y)\) is a smooth, bounded, function of x, y, and x' 
for any m. Then, if /3 > 0, 



<C( {ix',y)f^\^{x-x'mx,y)\-\{x',y)\)\\f{x',y)fdx') {x - x'f^\^{x' - xmx,y)\ - \{x' ,y)\)\dx' 



Since J{x — x')^'^\'^{x' — x){\{x,y)\ — \{x' , y)\)\dx' < C, where we allow the constant C to change from line 
to line, we have 



where for the last inequality we used that | / — x){\{x,y)\ — \{x' ,y)\)\)dx\ < C. A similar argument 
can be used when /3 < 0, using instead in the first step that for /3 < 0, {z)^ < C{x — .x')^'^((i', y))'^. 

A similar argument works for a derivative in the z/\z\ direction, and the argument can be iterated to get 



Lemma 6.9. Let (f>,(f>i G C~(M), with (1 - </>)(! - ^i) = 1 - and (j)i{y) = 1 if \y\ < yn + I- If 

- 0(y)).g e {z)-PL\W-) for \a\ < l< p - 1/2, <^(j/)ff e {zy^L^W'), x e Cr(Sr'), and * e 
Cb°°(R"-i) with supp* n : |C| = |A|/c±} = 0, then 

Px{zl\z\){i - <i>M)^{D.){^ - (A - iQfc^^r'g G {zf'^+'L\W^) 

for every e > and P e Diff^, with I < (3 - 1/2. 



Proof. We give the proof for ^ supported in y/\z\ > 0; the proof for x with support in y/\z\ < is similar. 

We use a cut off- function, \I>i e C^(R) with = 1 when |^| < A/c+ and supported in a small 

neighborhood of that region, so that supp(^'i(|^|)^'(^)) C {|^| < A/c+}. We write 





the lemma. 



□ 



(1 - My))x{z/\zmD.){A - (A - iO)%'r'g 




(39) 



The main contribution is 



(27r)i-"(l - My))x{z/\z\) j e-%(0*i(|^|)e-«(^^ + \i? - (A - 'Md^- 



(40) 



Here we may write, for y > yM and |^| < A/c+, 



{Dl + ICI^ - cf{\ - iQf)-'g{^, 0(2/) = e-'^V^V<-liP^(^) + ^^(^^ 



(41) 
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where £>^'5i G L2(R^"1; (y)-'^+^L'^{Ry)) for all k. Putting the first term of (|l|) into @, we obtain 

Note that if we apply + i\/c+, then the integrand vanishes on the critical set of the phase function. 
Since on the support of ^'^'i, g E H\ I a.n integer with I < (3— 1/2, we can integrate by parts to see that we 
have an element of 

For the tangential derivatives (in the z/\z\ directions), notice that if we have a derivative in a direction 
orthogonal to y, it commutes with {{D'^ + — Cq ^(A — iO)^)~^. That is, for example. 

By the decay properties of g and the regularity and decay properties of (1 — <t>)g, this gives an element of 
(z)i/2+eL2(K") after multiplication by x(l - </>i), if /3 > 3/2. 

After applying a derivative of the form t/g|- ~ 2;^^, as in the radial case the integrand vanishes on the 
critical set of the phase function, and so we can integrate by parts. 

This argument can be iterated up to / < /3 — 1/2. 

The second term of ( pl| ) gives an element of (z)^''+'+^L^(M"), where I is the order of the derivative P. 

On the support of (1 — ^'i(ICI)), ((-Dy + ICP) — (A — iO)^Cg ^)~^ is a smooth function of |^|, except near a 
finite number of points for which A^ is an eigenvalue of Cg (D^ + 1 ^ P ) • Since the eigenfunctions of this operator 
are exponentially decreasing in y and x is supported in?//|z|>(5>0 for some (5 > 0, these eigenfunctions 
do not contribute to the asymptotics here. Projecting off the eigenfunctions, we have 

(1 - ^,){Dl + - (A - zO)2cg-2)-ine : {y)-^L\^y) ^ {y)-^L\Ry) 

with bound C(|CP - C)-2. Therefore, 

P(l - <t>,{y))x{z/\z\) j ^i{\mDl + - (A - *0)2c-2)-\9(C, ^(yK e {z)-P+'L^{W^) 

where we used the fact that the inverse Fourier transform is an isomorphism on L^, and the regularity 
properties of £)"(1 — □ 

Lemma 6.10. Let f e (z)-°°L2(R"), (1 - (t>iy))f G 5(K") for some (j) e C^(R), and let u = (A - 
(A - iO)2)c~2)/. If^g e is m a neighborhood of \^\ = A/c± and xo e C;?°(§c); ^/len t/iere is a 

Wo = x± E/Jo* |2r^"'^+^^'j*o(£'x)M, wit/i e Diff^^(]R"), sMc/i </ia< 

(A - AVc^)wo = xo(2/|z|)V^*oPx)m + eo 

where eg G (z)-2./+5/2+ei2(]^n-)^ 

Proof. Let xo = X+ + X-i with x± supported in ±y > 0. We will outline the proof for xo = X+: the proof 
for Xo — X- is similar, and the functions can be added to get the general case. 

Recall that on the support of x+, V ^ ^j>J 1^1 jfy)- ^^'^ ^ ~ ^j=o ^^i- 
Let 

woo = (1 - My)M~'+'vjj^^^^x+'foiD,)u 
with G C^(]R), = and (piiy) = 1 for \y\ < yM- Then, by Lemma ^^6] 

(A - XVcDwoo = \z\-'vj{-^)x+'i'o{D^.)u + \z\-'-^P['fo{D.,)u + |z|-^+i*g(D,)(Fu + + 
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where P[ e Diff,^,(R") and G {z)~°^ L^{«"), so that (A - X^/cDwqo = x+V'^iDx)u + cqo, eoo e 

In the same manner, we choose wqi so that 

(A - XycDwoi = \z\-^-^vj+i{-^^)x+^oiDx)u - \z\--'~^P[^o{Dx)u + \z\--'-^ P^-^o{Dx)u + eoir 

with e Diff^(R") and eoir G {z)-^''+^/^+''L^{R"). Then (A - \^/4){woa + woi) = Vx+^^iD^) + eoi, 
eoi e (2:)^-'^'^/^+*L^(R"). This can be continued, with woj removing the terms in (z)^'^+^/^+'^^-'L^(M"), 
modulo terms in {z)--^-^/'^+''-i L'^{W'), up to j = J - 3. □ 

We shah also need the following lemma. 

Lemma 6.11. Ifw{x,y) e (zy/^+^L^iW), ^ e Cj^{W'^'^), andsnppw{^,y)nsupp^ = 0, then'^Vw{^,y) S 

Proof. Observe that 



^{OVw{^,y)^^{0 J ] e-'^-^V{x,y)e''''-^w{f^,y)dr^dx (42) 

= E*(^)/ J'^j{tvm-Vj)-'e-'-<^-^^D,^V{x,y)w{v,y)d7jdx (43) 
i 

where ^'j is a partition of unity with ^ rjj on supp \1/ j . We may repeat this integration by parts as many 
times as desired. Since \D"V{x,y)\ < Ca(z)^'^^l"l , the lemma follows. □ 



Proposition 6.3. Suppose x C^iS"^-^), * € C^{R"-^) /las supp *n{C : |C| = |A|/c±} = 0, supp^TljC : 
1^1 = '^j(A)} = 0, J = 1,2,...,T(A) , and tf ix,y)/\ix,y)\ € - x), y > 0, then ±Aa;/(c+|z|) ^ 

suppVE". Moreover, suppose if {x,y)/\{x,y)\ e supp{l — x), y < 0, then ±Aa;/(c_|z|) ^ supp^*, and f G 
{z)-°°L^{W-'^), (1 - 0(?y))/ e 5(R") /or some e C^'(M). Then 

x(n)*(^.)(A - (A - zO)Vc^)- V = e-^l^l/^"x(^)kr("-i)/^(ao(^) + Od^r^)) 

1^1 fI fI 

Proof Recall that if w = (A - (A - iOf/c^)''^f, then m = (A - (A - iO) Vcg)"i (Fu + /). If ^'(0^^^^/(C, 2/) 
were in C°°(R""^ {y)-°° L'^(R)), then using 

(A - (A - zO)2co-2)-ig = (2vr)i-" 1 ((|Cp + - (A - *0)2co-2)-ig(e, •)) (2/)^^ (44) 

and stationary phase, we would be done. However, it is not clear that Vu should be in C°°(]R"^^; {y)^°° L^{R)). 
We shall show that x^(^x)w can be written as a sum of two terms: one vanishing faster than u at infinity, 
and another of the form (A - (A - iO)Vco)"^5fc, where *(05fc(f, y) S C'=(R"~^ {y)-''L'^{M.)), where we can 
make k as large as desired. Then (^) and stationary phase will finish the proof. 

To do this, we follow an iterative procedure. The first step has been done in Lemma ^T^. We will 
iteratively construct functions wi which have the property that (A — A^/cq)(w — w;) improves with increasing 
I in an appropriate sense. 

Let ^I'o, *2, •■• e C^(M"~^) be such that, for all i, ^'i = I on the support of ^i+i^i = ^i+i, and 
^'i satisfies the support requirements placed on in the statement of the Proposition. Let Xo^ Xij X2) •■• S 
C^(§"~^) be such that xoX = X and Xi+iXi — Xi^ i = 0,1,2.... Let Wq be the function constructed in 
Lemma ^.ICj for this 4*0 and xo • Using the notation of that lemma, let 

to = (1 - Xo)V-^o{Dx)u + f-eo + V{l- ■^o{D,))u. 
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Note that the first three terms are in (z)^^'^+^/^+'^i^(R"), where we use the support properties of xo and 
along with an integration by parts argument as in the proof of Lemma to obtain the result for the 
first term. Additionally, the support of the Fourier transform in the x variables of (1 — 'i'o{Dx))u is disjoint 
from the support of ^'i. Then, using Lemmas |6.9| and |6.11 , we have Pxi^i{Dx){A — (A — iO)^/c§)~^to € 
(zy^'^+'L'^iW), for aU P e DifrJ.(M"), I < 2J-4. Since u^wa + (A - (A - iO) Vc^)"4o, which can be seen 
by a modification of the proof of Proposition |6.l| and the uniqueness result (Proposition [4.2[ ), this in turn 
means that P'i'i{Dx)wo G {z)--'+^^^+'L^{W') for aU P e DiffJ,, Z < 2 J - 4, using Lemma |64 
Given Wqi^o above, we now itcratively construct wi for / > 1 such that 

(A - XycDwi = Vxi'fiiDx)u + e, = Vxi'fiiDx)wi^i + yx;*;(^.)(A - (A - iO)Vcg)~4/-i + e,, 

where e/ G (z)--'+^/^+''-"^''-'^'> L'^iR"") and U is defined by 

ti = V ((1 - xi)'i'i{Dx) + (1 - ^^liDx))) wi-i 

+ V{{1- ^liDx)) + (1 - Xi)'fi{Dx)) (A - (A ~ iQf/cl)~hi^i -ei + f. 

Then 

u = wi + iA-{X-iOf/cl)-Hi. 
Moreover, U = t\+t'{, where y) e C°°(M^'"\ (j/)-°°L2(M)), t'l e (z)-('+i)('^-2)--'+i/2+ei2(-j^«)^ 

a.nA P'^i+i{Dx)wi e {z)--'+^/^+^L^{W^) for aU P e Diff;"(]R"), m < (/+2)(J-2). Additionally, suppw,(z) C 
supp xi{z/\z\). 

Supposing that wi^i,ti^i are as above, we show how to construct wi. Since 

Pxi^iiDx)wi-, e (z)-^+3/2+e^2(jg„) Pxi^iiDxKA~{X~iOf/co)-hi^i G {zy/^+'L^{W') 

for all P e Diff j^'^^'"'"'^^-' , we can, just as in Lemma |6.10 , find wi = ^^'^q'^'""'^^'^^ wij so that 

(A - X^cDw, = yx;*i(^.)«^/-i + Fx/«'/(^.)(A - (A ~ iO)Vc2)-it,_i + ei 
with e, e (z)-"'+i/2+e-(/+i)(J-2)^2(]^n)^ Let t[ = 1/(1 - ■^iiDx))iwi + (A - (A - i0)2)-it,_i)); then 
^z+i(OiJ(C,y) e C°°(R"~\ (y)-°°L2(R)) by Lemma|oil. Since supp (z) C supp x/-i(^/|^|), by Lemma 
6^, y(l - x/)^'/(^:.)w;/-i e (z)-°°L2(M"), and 

y(l - x/)*;(i?x)(A - (A - zO)2co e (^)-('+i)(./-2)-i/2-./+.^2(j^„) 



using the support properties of 5*; and xi and an integration by parts argument as in Lemma 6.9. Thus 

t'{ = ti-t[e (z)-('+l)(J-2)-./+l/2+e^2(lgn)^ 

Note that 

wij = \z\-^-'+^Pij{l - <l>{y))^,{Dx)[wi^i + (A - (A - tOf/clyHi^i] = {zr^-'+^Pijil - ^{y))^i(Dx)u, 

(45) 

with Pij e Diff^(K"), (/) e (M). Since (A - AVc§)(m - wi) = U, we have 

u = w, + (A-(A-iO)Vc2)-4z. (46) 

Note that P-fi+i(Dx){A-{X-iO)^/4)-Hi e {zy/'^+'L^{«") for P G Difr('+2^('^"^) (R"). Using M), (|4|), 
and Lemma U, this in turn means that P*;+i(P>:r)w; e (z)3/2+^-"'L2(]^n) f^j. p ^ Difr['+^^('^^^(R"). 
Thus, for any I > 1, wi and ti can be constructed to have the desired properties. 

To prove the proposition, we use (|6|). Since U ^ t[ + t'{ with ^i+i{^)i[(^,y) e C°°(R^'"\ {y)-°^L^{M)), 

t'l e (2)-(i+l)(J-2)-,/+l/2+e^2(R„)^ j^j^yg ^'(0^/(C,y) e iJ^(R"-l;(2/)-('+l)('^-2) + l/2+e-,/+._^2(]g)) fo^. 

s < (Z + 1)(J — 2) + J — 1/2 — e. Then, using equation ( ^ ) and stationary phase, we see that 
X*(i^.)(A - (A - iQf/cl)-hi = xe-^'l^l/l^r^"-'^/'(«o(^) + 0(Nr^)) 
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with ao e C('+i)(-''-2)+.-'-"-4(-§j ^hen I is sufficiently large. 

To finish, then, we only need show that '9{Dx)wi is of order |z|~'^"+^^/^. But we recall that 

|z|(^ + zA/c)d:/|,|^, e (z)-^+3/2+^l2(M"). 

This in turn means that '^{Dx)wi has the same properties, and that = C'(|z|^("^^^/^^'^+^+'^). □ 

6.4. The structure of the scattering matrix. Combining the construction of the approximate Poisson 
operator of Section^ and Theorem 1.3, we have now proved Theorem 1.1. For completeness, we remark that 
the construction of the approximation of P{\)Iij, 1 < j < T{X) in Section 5.2, Corollary 6.1, and Theorem 
1.3 prove 



Proposition 6.4. Let c,cq satisfy the general conditions of Section ^ and either hypothesis (HI) or (H2). 
Let A[X) = {Aij{\)), <i,j < T{X). Then, for j > 0, Ajj{X) is a Fourier integral operator associated with 
the antipodal mapping on S"^^, and Aij{\) is a smoothing operator when i ^ j- 

7. The Inverse Problem 



We recall our central inverse result. Theorem 1.2: 

Theorem . Suppose c and cq satisfy the general assumptions of Section ^ as well as either hypothesis ( HI ) 
or (H2), and n > 3. Then, if c+ = c_, the asymptotic expansion at infinity of c — cq is uniquely determined 
by Co and the transmitted singularities of the main part of the scattering matrix at fixed nonzero energy. If 
c+ < c_, then the asymptotic expansion is uniquely determined by cq and the reflected singularities of the 
main part of the scattering matrix at fixed nonzero energy. 

In proving the results for the inverse problem, we use the techniques of ^ and much of their language. 
We recall the arguments from these papers, noting the adjustments that must be made for the stratified 
case. 



Theorem 1.2 follows from the following theorem, which is somewhat stronger. 

Theorem 7.1. Suppose n > 3, ci and C2 satisfy the general requirements of Section^, and either (HI) or 
(H2), for the same cq. Let S'i(A), S'2(A) be the corresponding scattering matrices for some A £ M\{0}. //, for 
c+ = C-, the transmitted part of the main part of Si{X) — S2{X} is of order —I, then Ci{z) — C2{z) = 0(|z|^'^^). 
If for c+ < c- the reflected part of >5'i(A) — S'2(A) is of order ^l, then ci{z) — C2{z) — 0{\z\~''~^) . 

Proof. Let S*^) ^.(M") be the set of functions which, for \z\ > 1, are of the form J2j>o 1^1*^ (^/l^l); where 
a,(z/|z|)eC,-(§"-i\{(U,0)}). _ 

Suppose that ci — C2 — W £ S'^[^(K") and that the scattering matrices associated to cfA and c^A have 
the same transmitted (if = c_) or reflected (if c+ ^ c_) singular parts, to order I > k. Then we shall 
show that actually ci - C2 G 5'^^"^(M"), and thus by induction ci - C2 G S'^'7^(M"). 

If ci - C2 e S'cTt(I^")' then X'^{c^^ - c^^) = A^q ^^-^^ci + C2)(c2 - ci) = U, with 

j>k j>k 

and U = 0{\z\-''), where C/_j- ± E C^(El-^). Let 

U-j.+ {u;,u;n) ifa;„>0 

U-j-{lJ, UJn) if tiJn < 

Note then that the first fc — 2 terms in the construction of the Poisson operator carried out in Section |^ are 
the same, and the difference comes in the k — 1st term. 

Although many of the underlying techniques are the same, we shall treat the cases c+ = c_ and c+ ^ c_ 
serially. 
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W-j{u,LL>n) 



We begin with the case = c_. In the construction of the Poisson operator, the transmitted parts (that 
is, in the lower hemisphere) of the k — 1st terms differ by 



2Ac4 



''. T+{u;n,X) / W-k{s',0;u;){smsT-^ds', (47) 



almost as in (2.3) of |18). Here T_|_(A,tJ„) is the transmission coefficient determined by equations (|14|)-(16). 
We remark that in case c+ — C-, T+—T-. The transmission coefficient must be nonzero for A £ R, A 7^ 0, 
^ 0, < < 1. Therefore, as described in ||lj, Section 4], we can recover from the difference of the 
transmitted parts of the scattering matrices 

W-k{s,e]uj){sinsf-'^ds 



as long as a; = (oJ, uJn) with a;„ ^ 0. 

If the transmitted parts of the two scattering matrices are the same to order /c — 1, then 



Ik = / W-k{s,e-uj){smsY''ds = Q. (48) 

Since this is true for all u) with Un 7^ 0, we can differentiate with respect to the starting point twice, use 
sin^ s + cos^ s — 1, and find that Ik-2 = W-k{s,6;ijj){sm.sY~'^ds = as well (see [|9[ Section 5]). 
Therefore, if k is even, we reduce eventually to the case with fc — 2 = and if k is odd, to fc — 2 = 1. When 
k is even, differentiating one more time with respect to the starting point shows that W-k is even; for odd 
k two more differentiations show that W-k is odd. 
When k is even, we obtain 



W-k=Q 

'7 

for each closed geodesic 7 by joining together two half-geodesics. However, by ||l5|. Theorem 4.7], for n > 3 
the x-ray transform on S"~^ with domain restricted to smooth even functions is 1-1. Although W-k may 
have a jump discontinuity at uon = 0, it is smooth elsewhere. As in the proof of Corollary 4.19], by first 
taking a convolution with W-k we may assume that W-k is smooth and, applying the theorem, it is thus 0. 

If k is odd, we consider j^W-k, which is even. Since for each geodesic beginning at Zn = 0, sins is a 
constant multiple of Zn/\z\, we obtain j^W-k = for each geodesic 7, and again W-k = 0. 

When c+ ^ c_, we use the reflected singularities in the inverse problem. Recalling that W G S^^'l, we can 
recover from the reflected singularities, when tj„ > 0, 

R+{ujn,X) (^J^° W-k{s,e;iu){smsf-^ds + W-k{s' ,9-u;){sins')''-^ds'}j (49) 

The first integral is along a geodesic originating at u; and continuing to {{(p,0)} C S"^^; the second integral 
is along the reflection of the first geodesic when it meets (/>„ = and the path of integration ends at the 
point (— ZZJ, aj„). The variable s' is the distance from the point (uJ, — a;„). 

It is, however, more convenient to think of the sum ([l9[) as the single integral 



R+{uJn,X) / W-k+{s,e;iu){sms)''-^ds (50) 

where 

"^-'^^^^-{w-ki^,-^.) if<^„<0 (^1) 

and s is the distance from uj. It is fairly straightforward to see by symmetry that (^9|) and ( |50| ) are the same. 
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If we can show that ( pO| ) is sufficient for recovering W-k.+ {s,0;Lu){siTis)'^ ^ds for ah uj with w„ > 0, 
then the analysis used in the case c+ = c_ wiU show that W-k,+ is if the reflection coefficients agree to 
order —k. 

It suffices that i?+(a;„, A) is for at most an isolated set of a;„ with < < 1, for we can obtain 
the integrals for these isolated values of ujn by continuity. We recall that for < a;„ < ^1 — /ct , 
|i?+(cLi„, A)| = 1 ([^ Chapter 3]). Moreover, because co(y) — c± is compactly supported for ±y > 0, for 
fixed A e M, A) can be extended to a meromorphic function of Un in a neighborhood of < a;„ < 1, 

except near a;„ = ^1 — c^/c^, where it is a meromorphic function of (1 — c^/c^ + c^w^/c^)^/^. Therefore, 

its zeros are isolated, and we have shown that it is possible to recover W-k{iij) for w„ > 0. 

A very similar analysis works for the lower hemisphere, proving the theorem. □ 

We remark that this proof shows that if ci — C2 G S'ci(S"^^), then the main part of 5i(A) — 52(A) is of 
order — fc + 1. 

Corollary 7.1. Let c, Cq satisfy the general conditions of Section^ and either (HI) or (H2), and let n > 3. 
Then c+, c_, and the main part of the scattering matrix at nonzero fixed energy determine c modulo terms 
vanishing faster than the reciprocal of any polynomial at infinity. 

Proof. We need only show that cq can be recovered from c-|_ , c_ , and knowledge of the scattering matrix at 
fixed energy. The leading order singularity of the scattering matrix A{\) determines and is determined by 
i?±(A,a'„), r±(A,w„), A and c±, where T_|_ are defined by equations ( p^[]l6| ), and a similar definition 
gives r_, just as in one-dimensional scattering theory (see e.g. [||). 

Fix A e M \ {0}. We can think of equation (|lj) in the slightly more general form 

[Dl-\^{l/cl-l/c^_,)^e)^^Q, (52) 

a Schrodinger operator with potential — A^(1/cq — l/c+) which is either compactly supported (if c-|_ = c_) 
or "steplike" (if c+ < c_). We can define the reflection and transmission coefficients i?±(fc), T±{k), for ( |5^ ) 
as usual, as in (|T|-|l^, and R±{k) = i?±(A,tj„), f±{k) = T±{\un), when k ^ Acj„/c+. Moreover, R±, T± 
are meromorphic functions of k £ C if cj^ — c^, and if 0+ < c_, they are meromorphic functions on Z, the 
Riemann surface on which k and (A;^ — A^/c'^ + A^/c'i)^/^ are single- valued holomorphic functions. Therefore, 

knowing i?_|_(A, cj„), T+(A,a;„) for < cj„ < y^l/c^ — 1/ct determines ^+(fc), T+{k) on the whole plane (if 

c+ — C-) or Z {if c+ < c_). This in turn determines the eigenvalues of D^ — A^(1/cq — l/c'^) and the norming 
constants. These, together with c± and R+, determine the potential — A^(l/cg — l/c+) (e.g. ||Tl|, |^). □ 
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